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Abstract 

We investigate tiie interaction between the cosmological relic neutrinos, and primordial gravi- 
tational waves entering the horizon before the electroweak phase transition, corresponding to ob- 
servable frequencies today vq > 10^^ Hz. We give an analytic formula for the traceless transverse 
part of the anisotropic stress tensor, due to weakly interacting neutrinos, and derive an integro- 
differential equation describing the propagation of cosmological gravitational waves at these condi- 
tions. We find that this leads to a decrease of the wave intensity in the frequency region accessible 
to the LISA space interferometer, that is at the present the most promising way to obtain a direct 
detection of a cosmological gravitational wave. The absorbed intensity does not depend neither 
on the perturbation wavelength, nor on the details of neutrino interactions, and is affected only 
by the neutrino fraction f^. The transmitted intensity amounts to 88% for the standard value 
fu = 0.40523. An approximate formula for non-standard values of is given. 
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I. INTRODUCTION 



The presence in the Universe today of a stochastic background of gravitational waves 
(GWs) is a quite general prediction of several early cosmology scenarios. In fact, the pro- 
duction of gravitons is the outcome of many processes that could have occurred in the 
early phases of the cosmological evolution. Notable examples of this kind of processes 
include the arnplification of vacuum fluctuations in inflationary 1] and pre-big-bang cosmol- 
ogy scenarios P], phase transitions P], and finally the oscillation of cosmic strings loops P|. 
In most of these cases, the predicted spectrum of gravitational waves extends over a very 
large range of frequencies; for example, inflationary expansion produces a flat spectrum that 
spans more than 20 orders of magnitude in frequency, going from 10^^^ to 10^ Hz. 

The detection of such primordial gravitational waves, produced in the early Universe, 
would be a major breakthrough in cosmology and high energy physics. This is because 
gravitons decouple from the cosmological plasma at very early times, when the temperature 
of the Universe is of the order of the Planck energy. In this way, relic gravitational waves 
provide us a "snapshot" of the Universe near the Planck time, in a similar way as the cosmic 
microwave background radiation (CMBR) images the Universe at the time of recombination. 

The extremely low frequency (ELF) region (z/q < 10^^^ Hz) in the spectrum of primordial 
gravitational waves can be probed through the anisotropics of the CMBR. In particular, 
gravitational waves leave a distinct imprint in the so-called magnetic or B-modes of its 
polarization fleldj^, 0]. The amplitude of the primordial spectrum of gravitational waves 
is usually parameterized through the tensor-to-scalar ratio r, i.e., the ratio between the 
amplitudes of the initial spectra of the tensor and scalar perturbations in the metric. The 

Q n 

Planck satellite [7| , scheduled for launch in 2007, is expected to be sensitive'8" to r > 0.05. 
The lower limit corresponds to a density parameter Qgw{^) = {^/ Pc)dpGW / dlogu as faint 
as ~ 3 X 10~^^h^^ {h is the dimensionless Hubble constant) in the ELF range. Although this 
value looks incredibly small, it should be noted that, in order to produce such an amount in 
the framework of inflationary models (that at present time represent the most promising way 
to produce a signal in the region under consideration), a very early (starting at t ~ 10^^^ 
sec) inflation is required, and this possibility looks, from a theoretical point of view, quite 
unnatural. Future experiments are expected to enhance the sensitivity of two orders of 
magnitude P 
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On the other hand the planned large scale interferometric GW detectors, although de- 
signed with the aim to detect astrophysical signals, can possibly also detect signals of cos- 
mological origin 9|- They give complementary information with respect of the CMBR po- 
larization field since, even if their sensitivity is by no means comparable to the one than can 
be reached by CMBR polarization experiments, nevertheless they probe a different region in 
the frequency domain that would not be accessible to those ones. In particular the ground- 
based interferometers, such as the LIGoQ, VIRGO 1^, GEO600 12| and TAMASOoQ 
experiments, operate in the range 1 Hz < z/q < 10^ Hz, and are expected to be sensitive 
to VLGwh"^ > 10^^. Even more interesting is the LISA space interferometer [l^ . that will 
probably operate from 2013 to 2018. Not being hampered by the Earth seismic noise, it will 
probe the frequency region between 10~^ and 1 Hz and will in principle be able to detect 
^GV/h'^ ^ 10~^^ at z/q = 10~^ Hz. According to theoretical predictions, a large enough GW 
signal at this frequencies can be produced, with the appropriate choice of parameters, by a 
pre-big-bang accelerated expansion, by the oscillation of cosmic strings, or by the electroweak 
phase transition occurring at T = 300 GeV. 

In order to compare the theoretical predictions with the expected instrument sensitivities, 
one needs to evolve the GWs from the time of their production to the present. It is usually 
assumed that gravitons propagate in vacuum, i.e., they freely stream across the Universe. 
In this case, the only effect on a propagating GW is a change in frequency (corresponding 
to the usual redshift of the graviton energy caused by the expansion of the Universe), while 
the intensity of the wave remains the same. However, GWs are sourced by the anisotropic 
stress part of the energy-momentum tensor of matter, so that the vacuum approximation is 
well-motivated only when this can be neglected. It is already known that the anisotropic 
stress of free streaming neutrinos acts as an effective viscosity, absorbing gravitational waves 
in the ELF region, thus resulting in a damping of the B-modes of CMBrQ Q, H Q- 
The aim of this letter is to understand if something similar happens in other regions of 
the frequency domain. In particular, we consider GWs that enter the horizon before the 
electroweak phase transition (EWPT). This corresponds to an observable frequency today 
z/q > 10~^ Hz, i.e., to all waves possibly detectable by interferometers. In order to study 
this issue, we solve the Boltzmann equation for the distribution function of cosmological 
neutrinos and use its formal solution to obtain an integral expression for the anisotropic 
stress of neutrinos. We follow the approach outlined in Ref. [l^, but we have to take 
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into account the interactions of neutrinos with other background particles. The anisotropic 
stress is then inserted as an integral source term in the Einstein equation for the evolution 
of tensor metric perturbations. In the limit of very short neutrino mean free path, relevant 
for waves well below a frequency of 10^ Hz, this equation can be cast in purely differential 
form and a numerical solution be sought with standard methods. We find that the intensity 
of GWs is reduced to ~ 90% of its value in vacuum, its exact value depending only on one 
physical parameter, namely the density fraction of neutrinos. Neither the wave frequency 
nor the detail of neutrino interaction affect the value of the absorbed intensity, resulting in 
an universal behaviour in the frequency range considered. 

The paper is organized as follows. In section |n] basic statements are provided in order 
to formulate the dynamical problem; in particular the coupled Einstein-Boltzmann system 
is discussed in correspondence to the interaction between GW and neutrino plasma. More- 
over, a detailed analysis of the perturbative anisotropic stress tensor together with the optical 
depth is presented. In section llllj we develop the real calculation of the GW absorption 
by neutrinos, by virtue of an approximation scheme of the source term, based on the fast 
interaction rate of neutrinos. The time behaviour of the transmitted GW intensity is nu- 
merically determined, and its dependence on the neutrino content of the Universe outlined. 
In section IIV[ concluding remarks follow in order to stress the phenomenological relevance 
of our analysis, with respect to GW interferometric detectors. 

II. PROPAGATION OF COSMOLOGICAL GRAVITATIONAL WAVES 
THROUGH A NEUTRINO PLASMA 

A. Basic formalism 

We shall use, all throughout the paper, natural units in which c = h = = 1- 

Let us consider a gravitational wave, propagating on the background of a fiat Friedmann 

Universe. In synchronous gauge, the spatial components of the perturbed metric are written 

as: 

9ij = [Sij + hij] (1) 

while the other components are left unperturbed: goo = 1 and goi = 0. Here a{t) is the 
cosmological scale factor, and h^j is a small tensor perturbation representing a GW, that 
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evolves according tojl^: 

where TTj^ is the anisotropic stress part of the energy-momentum tensor T^^. It is defined 
through the relation vr^j = T^j — V6ij, where V is the background pressure. The space part 
of the energy momentum tensor can be computed by means of: 

T) = -^j f{x\ p„ t)^ dp, dp2 dp3. (3) 

where g is the determinant of the metric, p^ = dx'^/dX is the four-momentum (A is an affine 
parameter), and f{x\pj,t) is the distribution function (DF) of the particles in the fluid. 
The DF evolves according to the Boltzmann equation: 

L[f] = C[f] (4) 

where the L = df/dX is the Liouville operator, and C is the collision operator accounting 
for collisions between particles. Using the geodesic equation, the left-hand side (LHS) of 
the Boltzmann equation (the so called Liouville operator L) for / can be cast in the form 
(summation over repeated indices is understood): 



(5) 



dt p^ dxi 2p^ dxi dp. 
Once the collision term is also specified, eqs. (j21), © and (jH), together with the Friedmann 
equation describing the background evolution, are all that is needed, at least in principle, 
to follow the propagation of a GW. 



B. Anisotropic Stress Computation 

In this letter we shall consider GWs that enter the Hubble horizon well before neutrino 
decoupling (occurring at T ~ IMeV), and we shall suppose that the main contribution to 
the anisotropic stress of matter is given by neutrinos ^. We can write the DF of neutrinos 
as: 

f{x\p,,t) = f{x\p,M) + 5f{x\p,,t), (6) 

^ This approximation relies in the fact that, roughly speaking, the anisotropic stress is proportional to the 
mean free path of the particles in the fluid. In the framework of the standard model of particle physics, 
neutrinos are the most weakly interacting particles in the Universe, and then their anisotropic stress is 
larger with respect to other particle species. 
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so that at the initial time ti we have 6f = 0. We take for the DF at the initial time the 
usual Fermi-Dirac form ^: 

f{x\pj,h) = [exp {,/g^Jp-p]/T) + l]^^^^ , (7) 

where N counts the number of internal degrees of freedom. 

In respect of the case of GWs entering the horizon after neutrino decoupling (as in Ref. 
the situation we are going to deal with is complicated by the fact that the right-hand 
side (RHS) of the Boltzmann equation does not vanish, neutrinos being yet interacting with 
the cosmological plasma. The calculation of the exact collision term is usually not easy to 
deal with; however, a reasonable guess is given by assuming that C*[5/] = —p^Sf/tc, where 
tc is the mean time between collisions, and the p^ factor comes from writing the Boltzmann 
equation in terms of the conformal parameter A. 

The following relations will be useful in the following (we define p = ^JPiPi) : 

pO = ^ + 0(%), p^ = ^ + 0(M, (8) 
a a 

We Fourier transform the spatial dependence of the involved quantities. The comoving 
wavenumber ki is the Fourier conjugate variable to the comoving coordinate x*. Introducing 
the time variable u: 

, /■* dt' 

the Boltzmann equation finally reads, to first order in perturbed quantities: 

5f - ipikiSf - ^-PiPjPkh [hjk{.u) - /ijfc(O)] = ■ ^ (10) 

where the dot denotes derivative with respect to m, and hats denote unit vectors; /*^°^ is the 
zeroth-order part of /. Defining: 

J {u, ki, pj) = -ipih + (11a) 

aij[u, ki,pj) = —kipiPiPj (lib) 



^ Please note that this is not the usual splitting into a zeroth-order, equilibrium part plus a first-order part, 
since /(ti) contains itself a first-order part. However, Sf is certainly first order since the zeroth-order part 
of the DF O satisfies the coUisionless Boltzmann equation over the flat, unperturbed Robertson- Walker 
metric. 
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a formal solution for 6f satisfying the initial condition Sf{0) = can be written as (for 
shortness of notation we keep indicating only the dependence on u): 





pu 


nu 




6f{u) = exp 


- / J{u')du' 


X / exp 


r j{u")du" 




. Jo 


Jo 


Jo 



aij{u') hij{u') — hij{0) du. (12) 
Now the anisotropic stress tensor can be computed by means of equation (jHl); the result is: 



TTij = -Apu / < hijiu') + t{u') hijiu') - hij{0) 



X e 



t{u')—t{u) 



K{u — u)du, (13) 



where p^, is the mean (background) neutrino energy density, K is the kernel: 
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:i -x")e*"^ dx, 



(14) 



-1 



and r is the optical depth of the plasma with respect to the propagation of neutrinos: 



r(n) 



nat dt = — 
k 



ana du. 



(15) 



Here a is the relevant cross section for the reactions that keep neutrinos coupled to the other 
particles in the plasma, n is the density of target particles. We note that since na = tc, we 
have f = a/ktc. The dimensionless quantity r gives an estimate of the average number of 
interactions occurring between ti and t. 



C. Behaviour of the optical depth 

In order to proceed further, we need to calculate the optical depth t{u). In the following 
we shall consided GWs entering the horizon before the EWPT (T ~ 300 GeV). In order to 
better understand which region in frequency space we are investigating when we make this 
assumption, it is useful to find the relation between the gravitational wave frequency and 
the temperature of the Universe at the time when its horizon crossing occurs. This can be 
done very easily as follows. The temperature T^c at which a wave, having today a physical 
wavelength equal to Aq, entered the horizon is given by [2^: 

T.c(Ao)=63,;V2^,f (^y'eV, (16) 

where g^, and (7^,5 measure the contribution of relativistic species to the energy density 
and entropy density respectively. Before EWPT, all species in the standard models of 
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fundamental interactions are relativistic, resulting in = 106.75. Furthermore, in thermal 
equilibrium g^s = 9*- Then we find that the present day frequency z/q of the wave is: 

Quite interestingly, we find that all cosmological GWs that can be detected by laser in- 
terferometers, both spatial and ground-based, entered the horizon prior to the electroweak 
symmetry breaking. 

Neutrinos are kept in thermal equilibrium with the cosmological plasma through the 
weak interactions with other particles. At energies large as the ones we are considering, all 
the particles in the standard model of interactions are a possible target for neutrinos. The 
density n that appears in the expression for r is then 

C(3) 



n 



7r2 



-9*nT' (18) 



where ^(3) is the Riemann zeta function of order 3, and g^:^, in a similar way than g^, and 
(7*5, measures the contributions of relativistic species to the particle number density. At 
T > 300 GeV, g^N = 95.5. This yields n ~ 11.6 T^. 

For what concerns the cross section, we have that in the regime corresponding to average 
thermal energies much larger than the rest masses of the W^,Z^ weak bosons m^, rriz ~ 
100 GeV the temperature dependence of a is roughly given by: 

^ - (T » (19) 

where Gp = 1.1664 x 10~^ GeV~^ is the Fermi coupling constant. More detailed calculations 
can be performed to bring more precise formulas for the cross section. However, as we shall 
see in the following, the above expression us enough for our purposes. 

The last step in the computation of r is to express the comoving wavenumber k in terms 
of The. Considering the initial time ti to be very early, so that ti ~ 0, we have that in 
the radiation dominated era a{t) oc ^/t and a{u) oc u. Moreover, since /g dt'/a{t') is the 
comoving horizon size, we have that u = 1 corresponds to the time of horizon crossing. We 
choose to normalize the scale factor so that it is equal to unity at horizon crossing, then we 
take aiu) ^ Now we just need to cons.der the phys.ea. wavelength A as a funefon of 
temperature |20|] : 

.m.5.83xlO-(^).(^)".;V'GeV-. (20) 



Evaluating this expression at The and eliminating Aq by means of eq. (fT^. we get for the 
comoving wavenumber k: 

A; = 2.76- 10^1 (^^^Mpc-^ (21) 

Using this expression and eqs. (fTHj) and (fT^. together with fact that the quantity aTgl^^ is 
constant we find, for T ^ 100 GeV: 



III. ABSORPTION OF COSMOLOGICAL GRAVITATIONAL WAVES 

Now we can finally deal with the GW evolution equation. Defining xi'^) = hij{u)/hij{0), 
the equation reads: 



x + 2nx + x = -24un' x ^ K{u- u') x' + r{x' - 1 



e^'-^du. (23) 



Here Ti. is the "Hubble constant" a/ a = u^^, fi, is the ratio of neutrino density to the 
total density (in the standard model, fy = 0.40523), and r is given by eq. (j22I)- We use 
the short-hand notation that primed quantities are evaluated in u', and unprimed ones are 
evaluated in u. The above equation is an integro-differential equation for x- K is supplied 



by the initial conditions 



id: 



(24) 

m = 0. 



A numerical solution can be found by using standards method for integro-differential 
equations!^, 13]. It is however useful to consider two limiting cases. For waves enter- 
ing the horizon very early (i.e. The ^ 10^^ GeV, or equivalently z/q ^ 10^ Hz), the derivative 
of the optical depth is very small and then Eq. p3|) recover the free streaming form found 
in Ref. [l2|. It is found in that paper that this results in an absorption of 35% of the wave 
squared amplitude. On the opposite, for waves having The ^ 10^^ GeV, the derivative of the 
optical depth is very large. This makes the integration of the evolution equation not trivial, 
since the integral on the right hand side receives contribution only from a very small region 
where e"^'""^ significantly differs from zero. Technically speaking, the integral kernel is said 
be "quasi singular" |2ll. [2^ . However, we can exploit this characteristic to our advantage. 
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In fact, when the following disequality holds: 

<|r| (25) 

X - 1 

we can approximate e^''^' with a Dirac delta function 5{t — t'). The above condition can 
be checked a posteriori after obtaining a solution. However, apart from the mathematical 
details, it is clear that this approximation relies on the fact that the characteristic time 
scale for the variation of the gravitational field is much larger than the mean time between 
collisions. The integral source term on the RHS of eq. ()23|) depends, in general, on the 
whole past history of the system. Nevertheless, when the condition expressed in ()25|) holds, 
the very fast interaction rate destroys any information on the past history of the system, so 
that the integral is dominated by local (in time) contributions. 

This condition can, on the other hand, tested a priori using a test function in place of the 
actual solution xi.'u)- A reasonable guess is to suppose that the solution we are seeking does 
not differ too much from the vacuum solution Xvac? i-e. the solution in absence of anisotropic 
stress. The solution of eq.(PH|) with vanishing RHS is Xvac(w) = sm{u)/u. We use this as a 
test function. Inside the horizon (u > 1), where the source term is really effective^, we have 
that |xvac/(Xvac ^ 1)| < 1- The condition then simply reads f ^ 1. 

Seen in another way, the expression of the vacuum solution in terms of its Fourier com- 
ponents is: 

X.ac{u) = \ j ^e-'^'^dv (26) 

so that only harmonics with "frequency" |f | < 1 give a contribution to this solution. The 
condition that the mean collision time is much smaller than the time scale for the GW 
evolution, reads then f ^ v > 1 . Comparing with eq. (j22I), we see that this is the case 
for all k modes entering the horizon before the EWPT, up to modes entering as earlier as 
The ~ 10^^ GeV. 

Thus, using the substitution: 

e^'- - 5{r' -r) = (27) 



the evolution equation reads: 

2 . ^ 8/, 

-x + x = 

u bu 



X + -X + X = -^(X-1) (28) 



Well outside the horizon, Xvac — and Xvac — 1, so that the source term is very small. 
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where we have dropped a term proportional to x/'^'- 

This ordinary differential equation, describing the evolution of a gravitational wave in a 
tightly coupled plasma, is, together with the more general equation and the expression 
for the anisotropic stress, the main achievement of this letter. We note that, once the 
condition is fulfilled, the GW evolution does not depend on the details of the interac- 
tion, since the optical depth, and consequently the interaction cross section, do not appear 
anymore. The only physical parameter related to the background conditions appearing in 
eq. (UHl), is the fraction of energy density contributed by neutrinos. Moreover, using the u 
time parameterization, the equation is universal in the sense that it does not depend on the 
wavenumber k. 

A numerical solution to the above equation can be found quite easily. We plot this 
solution (let us it call the "matter solution") for the standard value = 0.40523, together 
with Xvac, in the left panel of Fig. [T] The solution follows quite closely the vacuum one. 
A closer inspection shows that, far after the horizon crossing, the squared amplitude of a 
wave propagating in the cosmological plasma gets damped. We use as a measure of the 
transmitted intensity the following quantity T: 

T ^ (29) 

\Avac/ 

where brackets denote time average over a period. In particular, we find that at m = 100, 
T 0.9, so that the intensity is reduced by a factor ~ 10% with respect to the vacuum case. 
We expect the damping not to become too much larger than this value, since the source 
term is multiplied by a factor 1/m^, and the gets rapidly suppressed after the entrance in 
the horizon. This is confirmed by examining the behaviour in time of T, shown the right 
panel of Fig. ^ We see that it approaches asimptotically a constant value, let us denote it 
with Too, equal to 0.88. 



Non-standard physics (for example, the presence of a lepton asymmetry) can result in 
a different value for f^y and consequently for the transmitted intensity 7^. Integrating 
the equation for different values of fu in the range [0, 1] we find that, with a very good 
approximation, 7^ depends linearly on f^. A fitting formula is given by: 

= 1 - 0.32/, + 0.05/J (30) 
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FIG. 1: {Left panel) Time evolution of the gravitational wave amplitude Solid line represents 

a GW propagating in neutrino matter. Dashed line represents a GW propagating in vacuum. {Right 
panel) Time evolution of the transmitted wave intensity T. 

IV. CONCLUDING REMARKS 



Our analysis outlined how gravitational waves produced in the early Universe are softly 
damped due to the anisotropic stress of cosmological background neutrinos. We have con- 
sidered gravitational waves entering the horizon before the electroweak phase transition, 
corresponding to observable frequencies today z/q > 10"'^ Hz. We find that the anisotropic 
stress is closely related to the optical depth r of the cosmological plasma with respect to 
the propagation of neutrinos, and to its first derivative with respect to time. In the limit 
of very small time derivative of the optical depth, relevant for frequencies uq ^ 10® Hz, 
the behaviour studied in Ref. Q| is recovered, resulting in an absorption of 35.6% of the 
wave intensity. This frequencies are however unaccessible to present-day interferometers. In 
the opposite limit of very large values of this derivative, relevant for observables frequencies 
uq ^ 10® Hz, we find that roughly 12% of the wave intensity is absorbed by the cosmological 
plasma, independently from its wavelength. This result does not depend on the details of 
the weak interactions between particles in the plasma, but depends only on the amount of 
neutrinos present in the Universe, as parameterised by the quantity f,^. 

The importance of our results relies in the fact that the damping affects GWs in the 
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frequency range where the LISA space interferometer and future, second generation ground- 
based interferometers can possibly detect a signal of cosmological origin. This effect is 
roughly of the same order of magnitude as the one affecting GWs detectable through the 
B-modes of CMB polarization. The damping is not so severe to make the detection of 
cosmological waves unfeasible by interferometers. However it should be taken into account 
when testing the theoretical predictions of early Universe scenarios against observations. 
Moreover, the dependence of on f^, can be exploited to measure the latter, and to 
constrain models of non-standard physics. This even more important in view of the fact 
that in this way we could measure the value of fy at very early times, while available 
constraints regard the neutrino fraction at the time of cosmological nucleosynthesis or at 
the time of matter-radiation decoupling. 
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